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Spectral components in the bright, long GRB 061007: 
properties of the photosphere and the nature of the outflow 
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ABSTRACT 

We present a time-resolved spectral analysis of the bright, long GRB 061007 (z = 
1.261) using Swift BAT and Suzaku WAM data. We find that the prompt emission 
of GRB 061007 can be equally well explained by a photospheric component together 
with a power law as by a Band function, and we explore the implications of the former 
model. The photospheric component, which we model with a multicolour blackbody, 
dominates the emission and has a very stable shape throughout the burst. This compo- 
nent provides a natural explanation for the hardness-intensity correlation seen within 
the burst and also allows us to estimate the bulk Lorentz factor and the radius of the 
photosphere. The power-law component dominates the fit at high energies and has 
a nearly constant slope of —1.5. We discuss the possibility that this component is of 
the same origin as the high-energy power laws recently observed in some Fermi LAT 
bursts. 
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1 INTRODUCTION 

The prompt emission of gamma-ray bursts (GRBs) is usually 
well modelled by a smoothly broken power law ( Band ct al. 
1 19931 ). peaking in the 100-1000 keV energy range. The phys- 
ical origin of this emission is still unclear, and the situ- 
ation has recently become even more puzzling, with the 
Fermi satellite revealing that at least some bright bursts 
have additional power-l aw components extending well into 
the GeV energy range ( Abdo et al] l2009l : lAckermann et al.] 
l2010l : lGranot" et al ] |201oi r For a long time synchrotron emis- 
sion has been considered the most promising candidate for 
explaining the Band component of the spectra. Although 
this model can explain many aspects of the emission, it also 
suffers from some significant problems. In particular, many 
bursts are observed to have harder low-energy spectra than 
predicted by standard synchrotron models (ICrider et al.l 
ll997l : |Preece et al.lll998l : lGhirlanda et al.ll2003l '). 

An alternative model that is able to account for 
the h ard low-energy s pectra i s photospheric emis- 
sion l|Meszaros fc ReesI l2000l : iMeszaros eT^ |2002| : 
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iRees fc Meszaroj l2005l : iPe'er et all 120061 '). Such a model 
also has the advantage of providing a natural explanation 
for the observed correlations between the peak energy 
and lu minosity within/be t ween b ursts |Golenetskii et akl 
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1 19831 : iBorgonovo fc Rvdd I2OOII : lYonetoku et~ 

Fits to observational data also support this picture, 
with blackbody models (often in combination with an 
additional non-thermal com ponent) providing excellent 
descriptions of many bu rsts (|Ghirlanda et al.l l2003l : I Rvdd 
I2OO4I : iRvde fc Pe'eill2009l ) . However, in spite of its successes, 
it seems clear that single, narrow blackbody components do 
not sig nificantly contribute to the prompt emission of most 
bursts (jGhirlanda et al.ll2007l : lBellmll2010l ). 

This does not necessarily rule out a photospheric origin 
for the emission. In fact, there are many reasons to expect 
the photospheric emission to deviate from a simple Planck 
function, and several authors have recently considered mod- 
els that give rise to photospheric emi s sion with a broad , 
Band-like spectrum (Pe'er ct al. 2006; 'Be loborodwl I2OIOI : 
iLazzati "fc^Bc gclman 2010; Mizuta ct al. 201^) . A broadened 
photospheric component is also to be expected from the 
simple fact that the observed spectrum is most likely com- 
posed of emission arising from different regions in space. As 
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a result, the angle dependence of the optical depth and the 
Doppler shift, as well as possible angle-dependent density 
and Lorentz -factor profiles, rnay all contribute to create the 
broadening (|Pe'er et al.ll2007l : IPe'er fc Rvdell2010l '). 

Here we analyse the bright, long burst GRB 061007 
with the aim of exploring models in which the emis- 
sion is dominated by a photospheric component. The 
prompt emission of GRB 061007 was caught by both 
Swift and Suzaku llOhno et ahlfioogl ) as well as Konus- Wind 
l|Golenetskii et al.ll2006l 'l. With an isotropic en ergy release 



of~ lO'^* erg in the 1-10 000 keV energy band (|Ohno et all 
l2009h this burst is among the most energetic GRBs ever ob- 
served, comparab le to the brightest bursts detected by the 
Fermi LAT (e.g. ICenko et aTJl2010l '). A very bright after- 
glow w as observed by Swift and several ground-b ased tele- 
scopes (|Mundell et al.1 l2007l : ISchadv et all 120071') . and the 
redshift was measured to be z = 1.261 ( Osip et al.l l200(j . 
[jakobsson ct al. 2006). 

Due to its brightness and the wide energy band of- 
fered by the combined Swift and Suzaku observations, 
GRB 061007 is an ideal candidate for detailed time- 
resolved studies of the gamma-ray emiss ion. The Suz a ku an d 
Swift data were previously analysed bv lOhno et all (|2009l ). 
who carried out a time-resolved spectral analysis on a 1-s 
time-scale. Modelling the spectra with a Band function, the 
authors found that the time-resolved spectra follow the same 
Speak — I/iso relation as the time-averaged spectra of other 
bursts (|Yonetoku et al.li2004 '). but also noted that the initial 
rising phases of the pulses may be outliers to this relation. 

Here we find that that an alternative interpretation in 
terms of photospheric emission is possible, and explore the 
consequences of this in terms of the evolution of the fireball 
properties. In addition to a thermal component, our best- 
fitting spectral model also requires the presence of a power 
law. We discuss the properties of this component in view of 
recent Fermi results. This paper is organized as follows. We 
describe the observations in section [2] and present the time- 
resolved analysis in section|31 We then explore the properties 
of the relativistic outflow as derived from our spectral fits 
in section [4] A discussion and our conclusions are presented 
in section O 




Figure 1. BAT + WAM light curves on a 1 s time-scale. The BAT 
count rates are per fully illuminated detector for an equivalent 



on-axis source. 



3 TIME-RESOLVED SPECTRAL ANALYSIS 

Fig.[T]shows the WAM and BAT light curves of GRB 061007 
on a 1 s time-scale. The light curve covers 70 seconds and 
three main pulses can clearly be identified in both energy 
bands. The time bins shown in the light curve are also used 
for the spectral analysis, excluding the first and last bins as 
well as the time between 19 and 28 seconds, due to the low 
count rates during these times. This leaves us with a total 
of 58 spectra with 1-s duration for the analysis. 

The spectra were fitted over the 14-150 keV (BAT) 
and 150-5000 keV (WAM) energy ranges, with the cross- 
normaliza tion between the t wo detectors left as a free param- 
eter (as in lOhno et al.|[2009l ). Fixing the cross-normalization 
at the value obtained from the time-averaged spectra only 
changes the results marginally, as discussed in section 13.41 
All fits were performed using xspec12, assuming a stan- 
dard cosmology with Q.x = 0.73, Q,m = 0.27 and Hq = 
71 km s~^ Mpc~^. Errors on model parameters and error 
bars in plots represent the one sigma confidence level for 
one interesting parameter, unless otherwise stated. 



3.1 Identification of a photospheric component 



2 OBSERVATIONS AND DATA REDUCTION 

GRB 061007 was observed with the Suzaku WAM and the 
Swift BAT detectors. The Tgo duration was measured to be 
59 s i n the 50-5000 keV WAM energy range (lYamaoka et al.l 
l2006l) and 75 s iii the 150-500 keV BAT energy range 
i Schadvet al]l2007D . The observations are described in de- 



tail in lOhno et alf l|2009l ), and for this work we also fol- 
low the data reduction procedure described in that paper. 
For the WAM analysis we thus use the WAM-3 detector, in 
which the burst was most strongly detected. For both the 
WAM and the BAT we extract science products on a 1 s 
time-scale, as this corresponds to the temporal resolution of 
the transient class WAM data. 



As a first step we investigate whether the spectra can be 
fitted with a model in which the peak of the emission is pro- 
vided by a single Planck function. To account for any non- 
thermal emission we also include a second component in our 
models. We consider two forms for this component, a power 
law and cutoff power law. For the majority of the spectra 
the blackbody-|-power-law model does not provide an ac- 
ceptable fit. Replacing the power law with a cutoff power 
law results in excellent fits, but we note that in this case 
the contribution from the blackbody is very small. This still 
holds if the energy of the blackbody is constrained in order 
to correspond to the peak of the spectra. We thus conclude 
that the prompt emission of this burst is not dominated by 
a single blackbody component. 

As mentioned in section [T] another, more realistic, pos- 
sibility is that the photospheric emission is broadened. We 
model this scenario with a multicolour blackbody, using the 
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XSPEC model diskpbb. This model was developed to describe 
emission from an accretion disc, but the parameters can also 
be interpreted in the general case of multicolour blackbody 
emission. The free parameters of the model are the shape 
of the blackbody (described by the parameter p) , the maxi- 
mum blackbody temperature (Tmax) and the normalization. 
More specifically, the parameter p describes the temperature 
profile of an accretion disc giving rise to multicolour black- 
body emission. Since this is irrelevant in the case of GRBs, 
we instead introduce the parameter q, which relates to the 
fitted parameter p as q = 4 — 2 /p. Using q we can express 
the relationship between the flux and the temperature of the 
single blackbody components (which together make up the 
multicolour blackbody) as 



F{T) = Fm 



\ -t max / 



(1) 



where Fmax is the flux of the Planck function at Tmax- As 
g — ^ oo we approach the case of a single blackbody. We 
note that the same parameter q is used t o des cribe the 
multicolour-blackbody model in lRvde et al.l (|2010f ). We also 
use ([T]) to define the average temperature of the multicolour 
blackbody 



1+1 

'q + 2 



(2) 



To allow for an additional non-thermal component we 
also include a power law in the model, N{E) = KE" , where 
K is the normalization at f keV in photons keV~^ cm~^ s'^ 
and s is the photon index. The resulting multicolour- 
blackbody + power-law (mBB+pl) model provides an excel- 
lent fit to all the spectra. Interestingly, we also find that the 
shape of the multicolour blackbody stays remarkably con- 
stant, pav = 0.70 ± 0.06 (corresponding to q — 1.1), based 
on the 47 spectra for which no problems occurred in the error 
calculations. We therefore re-fit the spectra with p fixed at 
0.70 in order to obtain better constraints on the remaining 
model parameters. Fig. [2] shows the resulting values of 
compared to those obtained from fits to a Band model (both 
models have 80 degrees of freedom). Clearly it is not pos- 
sible to distinguish between the two models on a statistical 
basis. The cross-normalization (C) between the WAM and 
the BAT is also similar for the two models (C = 1.27 ± 0.26 
for Band and C = f .f7 ± 0.19 for mBB+pl), showing that 
this parameter does not significantly affect the results (see 
also the discussion in section r3.4[) . Examples of fits to the 
mBB-|-pl model are shown in Fig. [S] 

The evolution of kT with time is shown in the top left 
panel of Fig. 2] A comparison with the light curve in Fig. [T] 
reveals that kT correlates with the fiux. We will explore this 
relation in more detail in section [531 below. In the lower left 
panel of the same figure we plot the dimensionless quantity 
{FBB/aT^y^'', which is a measure of the effective size of the 
emitting region. This parameter is seen to increase during 
the rising phase of each pulse and then only vary mildly 
during the decaying phase. In particular, the variations in 
this parameter are much smaller than the variations of the 
light curve (Fig[T]), supporting the interpretation that this 
burst is dominated by a photospheric component. 

In Fig. 3] we also plot the ratio of the blackbody fiux to 
the total fiux and the power-law photon index. We see that 
the contribution from the blackbody component is fairly sta- 
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Figure 2. Comparison of for the Band and multicolour- 
blackbody models. Both models have five free parameters. 
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Figure 3. Sample fits to the mBB+pl model. Plots in the left 
column are for the spectrum at time 16 s (at low fiux) and the 
right column are for the spectrum at time 41 s (at high fiux). 
The top row shows the count spectra and the residuals of the 
fit, with BAT data shown in black and WAM data shown in red. 
The bottom row shows the deconvolved components of the best- 
fit models. The solid, black line is the total model, the dashed, 
red line is the multicolour blackbody and the dotted, blue line is 
the power law. 



ble at around 75 per cent of the flux, although less well 
constrained during the first pulse. Due to the low flux of 
the power-law component in the fitted energy interval, the 
photon index could not be constrained in all the spectra. 
Considering only those 31 spectra for which it could be con- 
strained, we find that a fit to a constant gives Sav = —1.50 
with x^/d-o-f. = 30.8/30. 
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Figure 4. Results from the mBB+pl fits as a function of time. The shape of the blackbody was kept fixed at p = 0.70. Top left: evolution 
of the temperature of the multicolour blackbody. Bottom left: the effective size of the emission region. Top right: contribution of the 
blackbody to the total flux. Bottom right: power-law photon index. Data points without error bars indicate that the error calculation 
did not converge. 



3.2 Significance and properties of the power law 

It is interesting to consider the power-law component of 
our model in view of recent Fermi observations, which have 
revealed additional power-law compone nts at high ener- 
gies in a number of b right GRBs (e.g. lAbdo et al.l l2009l : 
lAckermann et al]|2010l ). We especially note that our pho- 
ton index values are clustered around —1.5, which has also 
been observed in several Fermi bursts ijGuiriec et al.|[2010al l 
and which is indicative of synchrotron or inverse Compton 
emission, as discussed in Section [5] However, we have also 
noted that the contribution from the power law to the total 
model is small and that the photon index cannot be con- 
strained in many of the time intervals, which prompts the 
question of whether the power law is required in the fits. 

To test this we fit all spectra with a pure multicolour- 
blackbody model (p still fixed at 0.70) and use an F-test to 
compare these fits with the mBB-|-pl ones. The results are 
plotted in Fig. [5] We see that the significance of the power 
law is well above 90 per cent in the majority of the spectra, 
but that the significance is lower at low fluxes. As expected, 
the photon index is also poorly constrained in these low- 
flux intervals (cf. Fig. U)). In order to make sure that the 
high significance of the power law in most of the spectra 
is not due to the fact that p was kept fixed, we also carry 
out the same test after fitting the spectra with a multicolour 
blackbody with p as a free parameter. In this case the overall 
significance of the power law is somewhat weaker, but it is 
still required at very high significance in the majority of 
spectra after the initial, weak pulse. 

Based on our spectral fits we also derive fiux hght curves 
for both the multicolour blackbody and the power-law com- 
ponent. The correlation between these hght curves was in- 
vestigated by calculati ng the discrete cro s s corr elation func- 
tion (DCCF) following lEd"elson fc Kroli^ l| 19881 ). The result- 
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Figure 5. Significance of the power law according to an F-test. 
Filled circles indicate spectra for which the power law is required 
above 99 per cent confidence. The dashed line is the light curve 
of the burst. 

ing DCCF is shown in Fig|6] The relative lag and the correla- 
tion strength were computed by fitting a Caussian profile to 
the correlation peak in the DCCF, and uncertainties in these 
parameters were esti mated using a Monte Carlo method (see 
IPeterson et aLlflQQSl ). The result of this analysis is a corre- 
lation strength of 0.44 +/-0.14 (corrected for the effect of 
measurement noise) and a lag which is consistent with zero, 
-0.8 +/-1.5 sec, where negative lag means that the power 
law leads the blackbody component. 

3.3 Spectral evolution during individual pulses 
lOhno et al.l l|2009l ) has already considered the iJpcak — ^iso 
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Figure 6. The cross correlation function between tlie blackbody 
and power-law components. Negative lag means that the power 
law flux leads that of the blackbody. 
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Figure 8. The photospheric radius versus the saturation radius 
of the fireball. Typical one-sigma error bars are shown in the lower 
right corner. 



evolution of this burst as derived from fits to a Band 
function. They found that the time-resolved spectra satisfy 
the i?Deak — Liao relation defined by time-averaged bursts 
ijYonetoku et al.l[2004 'l. but that a number of points in the 
initial rising phases of the three pulses are outliers to this 
relation. Here we instead consider the spectral evolution in 
terms of the flux and temperature of the multicolour black- 
body, which provides a natural explanation for the relation 
between peak energy and luminosity. 

The temporal evolution of each pulse in the Fbb — kT 
representation is shown in Fig. [T] The pulses exhibit a track- 
ing behaviour, with the rising phase of each pulse mov- 
ing along a lower track in the diagram. In addition, each 
pulse starts rising slightly higher up than the previous one. 
To quantify the temporal evolution within the pulses we 
fit straight lines to the logaritmized data in the rising and 
decaying phases of each pulse, using a fitting routine that 
accounts for the errors in both coordinates. We find that 
the correlation indices of the decaying phases are clustered 
arou nd 2.2-2.3, which is typical for decaying ph ases of GRBs 
(e.g. iBorgonovo fc Rvdell200]] . ILu et al.ll2010l ). The rising 
phases have steeper indices (values around 3) but are poorly 
constrained due to the small number of data points. 

As the three main pulses are clearly separated we can 
also carry out the experiment of treating them as individual 
bursts, and fit their time-averaged spectra with the mBB-|-pl 
model. The results of the three time-averaged fits fall within 
the range of results from the 1-s fits, as shown in Fig. [T] 
This is to be expected in the picture where the the various 
hardness-intensity correlations between the time-averaged 
spectra of different bursts are driven by the behaviour tha t 
is seen within individual bursts (see e.g lFirmani et al.ll2009h . 

3.4 Impact of the cross-normalization between 
the WAM and the BAT 

In the time-resolved fits described above we have allowed the 
cross-normalization factor (C) between the WAM and the 
BAT to be a free parameter. C was also l eft free in the prev i- 
ous analysis of this burst carried out bv lOhno et al.l l|2009l ). 
Since the Suzaku satellite kept the same attitude during 



the burst the C-values obtained from the time-resolved fits 
should in principle be consistent with the C-value from a 
time-averaged fit. However, we find that fixing C at the 
time-averaged value results in significantly worse in 
intervals with high fiuxes, regardless of the type of model 
being fitted. This is most likely due to systematic errors in 
the WAM detector, which increase with the brightness of 
the burst. 

As an alternative to keeping C free we therefore per- 
formed all the time-resolved fits with C fixed at the time- 
averaged value, but with 3 per cent systematic errors added 
to the WAM energy channels. The quality of the resulting 
fits is very similar to the fits with C free. The only systematic 
change in terms of best-fitting parameters is that p is slightly 
higher, pav = 0.73 ± 0.08 compared to pav = 0.70 ± 0.06 for 
C free, which translates into a slightly narrower multicolour 
blackbody. 



4 PROPERTIES OF THE FIREBALL 

Having identified a photospheric component in the spec- 
tra, we can now use the measured temperature and fiux 
of this component to determine the properties of the rela- 
tivistic outfl ow. Assuming the non-dissipative flreball model 
described in lPe'er et al. |(|2003) we calculate the bulk Lorentz 
factor of the flow (F), the radius of the photosphere (rph), 
the size at the base of the fiow (ro) and the saturation ra- 
dius (r-s). In order to determine all these parameters we must 
have > rs, wh ich we confirm following the procedure in 
IPe'er et al.1 (|2007l ) (see also Fig. H]). Fig. [9] shows aU the pa- 
rameters of the outflow as a function of time. 

r is seen to vary between about 200 and 600 (C^)^''*, 
where ^ is a geometrical factor of order unity and Y is the 
ratio between the total fireball energy and the energy emit- 
ted in gamma rays. The highest values of F are seen dur- 
ing the second pulse. The other three parameters plotted in 
Fig. |9] {rph, ro and Vs) all show a trend of increasing dur- 
ing the rising phase of each pulse and then remaining nearly 
constant during the decaying phase. The average values of 
these parameters are rph = 2.5 ± 0.8 x lO" C'^^^ Y^^'^ cm. 
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Time (s) 

Figure 7. Evolution of tlie flux and temperature of the multicolour blackbody during individual pulses. The first pulse is tracked in 
red in the top left panel, the second pulse is tracked in blue in the top right panel and the third pulse is tracked in green in the middle 
left panel. The middle right panel shows the results of time-averaged fits to the three pulses as a red square, a blue filled circle and a 
green triangle, respectively. The bottom panel shows the light curve with the three pulses identified. The initial rising phase of each pulse 
clearly follows a lower track than the rest of the pulse. In addition, we see that each pulse starts rising on a slightly higher track than 
the previous one. 
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Figure 9. Properties of the outflow, calculated from the best-fitting parameters of the mBB+pl model. ^ is a geometrical factor of order 
unity and Y is the ratio between the total fireball energy and the energy emitted in gamma rays. F is the bulk Lorenz factor, ro is the 
size at the base of the flow, rpji is the radius of the photosphere and is the saturation radius. 



ro = 3.5 ± 1.7 X 10* Y 



cm and = 1.3 ± 0.6 x 
10" ^-15/4 y-5/4 should be noted that there is some 

uncertainty in the last two parameters as they depend on the 
assumed fireball dynamics during the acceleration phase (see 
iPe'er et al.ll2007l for details). 

There is also a systematic uncertainty in our estimates 
of all these parameters due to the fact that they all de- 
pend on the total gamma-ray flux. For the results presented 
in Fig. [9l we have simply used the flux in the observed 
BAT-I-WAM energy band, but the flux will of course be 
higher if the power-law component extends to high ener- 
gies. The effect of a higher total gamma-ray flux would be 
higher values of F and rph and lower values of ro and Vs, 
with the change being significantly larger for the latter two 
parameters. 

In order to estimate to what energy the power law is 
likely to extend we consider the value of ro. As this is the 
initial size of the fireball we would expect a value of about 
10^ cm, i.e. the Schwarzschild radius of a solar-mass black 
hole. We therefore determine where the power laws from all 
our spectral fits would need to break in order to be consis- 
tent with this value. Assuming a break to a slope of —2 we 
find that the majority of the spectra should have a break 
in the 10 MeV - 1 GeV band, while a sharp cutoff would 
give breaks in the 100 MeV - 100 GeV range. If we also use 
the inferred power-law breaks to calculate new values of the 
other parameters, we find that F and rph are up to a factor 
of two higher than in Fig. |9] while Vs is between 10 and 50 
times smaller. 



5 DISCUSSION AND CONCLUSIONS 

After the launch of Fermi it has become increasingly ap- 
parent that GRB spectra are composed of multiple com- 



ponents. While some spectra are still very well described 
by the canonical Band function, others require an addi- 
tional power-law component extendin g to high energie s 
and/or a photospheric component (e.g. IZhang et allbOld ). 
The power law has b een detected with high significance 
in a h andful of burs t s (|Abdo et al. |l2009l: lAckermann et al.l 
l20ld : iGuiriec et al] l2010al : lAckermann et al] l201ll) but 
there is marginal evidence in many more i Granot et al.l 
I2OIOI ). A photospheric component is clearly seen in both 
GRB 090902B (together with a power-law component, 
iRvde et alj|2010l') and GRB 1 00724B (together with a Band 
function. iGuiriec et al.ll2010bt ). In the case of GRB 090902B 
the photospheric component is best modelled by a multi- 
colour blackbody. Given this evidence it is natural to ask 
whether such composite models can explain the prompt 
spectra observed with other satellites than Fermi . We have 
found that this is indeed the case for the Swift + Suzaku data 
of GRB 061007, which can be equally well fitted by a mul- 
ticolour blackbody together with a power law as by a Band 
function. Below we discuss the properties and possible origin 
of the photospheric and power-law components. 



5.1 The photospheric component 

We find that the multicolour blackbody dominates the 
prompt emission of GRB 061007, providing about 75 per 
cent of the total fiux throughout the burst. Since the 
luminosity and temperature of a blackbody are related 
this model also provides a natural explanation for the 
hardness- intensity correlation within the burst. This is 
a major advantage compared to models where the bulk 
of the emission is interpreted as synchrotron emission, 
where these correlations are much h arder to explain (e.g. 
iRamirez-Ruiz fc Llovd-Ronningir2002l '). 

Another interesting result from our fits is that the shape 
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of the multicolour blackbody is consistent with staying con- 
stant throughout the burst. The same shape is also found 
when fitting the time-averaged spectra of the three main 
pulses. This result disfavours a scenario where the broad- 
ening is primarily due to rapid temperature fluctuations, as 
this should give rise to a narrowing of the multicolour black- 
body whe n shorter time intervals are consider ed (as in GRB 
090902B, iRvde et al.|[2oTol : IZhang et al.ll2010l ). 

Several other mechanisms that may create broadened 
photospheric spectra have recently been discussed in the lit- 
erature. These include sub-photosphe ric heating due to c ol- 
lisions between neutrons and pr otons llBelo borodov''201(t) or 
shocks or magnetic dissipation (jLazzati fc Bcgclman 201(|). 
These models result in high-energy slopes that match those 
of typical Band spectra, but fail to explain the observed 
range of low-energy slopes. In order to explain the low- 
energy part of the multicolour blackbody in GRB 061007 
(which can be approximated by a power law with photon 
index a ~ —0.8) a different explanation is clearly needed. 
The most important effects that could contribute to a softer 
low-energy slope can be summarized as follows: 

(i) The observed temperature depends on the Doppler 
boost, which in turn depends on the angle to the line-of- 
sight (Tob = DT' = T' /T{l-l3cose), where T' is the outflow 
temperature measured in the comoving frame). Integrat- 
ing the emitting surface over ang le will therefore prod uce a 
multicolour blackbody spectrum. iPe'er fc Rvd3 (|2010l ) have 
shown that for late times, high latitude effects produce a 
spectrum with a = — 1. 

(ii) The assumption that all photon s originate f rom r(r = 
1) is not always a good approximation l|Pe'eij|200^ '). Depend- 
ing on how the comoving outflow-density scales with radius, 
the probability density function for a photon to make its last 
scattering at radius r does not have to be sharply peaked 
around r(r = 1). If the local outflow temperature or Lorentz 
factor varies on distance scales similar to the width of the 
probability-density profile, an integration over radius is nec- 
essary to obtain the observed spectrum. 

(iii) The outflow properties around the photosphere are 
expected to vary on a tim e-scale shorter than th e most 
highly time-resolved spectra (|Rees fc Meszaro3l2005f ). Thus, 
observed spectra are likely the result of an integration over 
time, which can produce a soft low-energy slope for rea- 
sonable temp oral scalings of the e mitting surface and the 
temperature (|Blinnikov et al.lll999l '). As discussed above, ob- 
servations indicate that this is not an important effect for 
GRB 061007. 

While all of these effects have the potential of softening 
the low-energy slope from purely thermal emission, simu- 
lations are needed in order to determi ne their relative im - 
portance. Indeed, recent simulations bv lMizuta et al.l (|2010l ) 
show thermal emission from a GRB jet with a — —0.5, much 
softer than a blackbody spectrum with a — 1. 

5.2 Properties of the outflow 

Having identifled a photospheric component in the spectrum 
we were also able to derive the properties of the outflow, in- 
cluding the Lorentz factor, the radius of the photosphere 
('"ph), the initial radius of the outflow (ro) and the satura- 
tion radius (rg). As already discussed, the latter two param- 



eters depend on the assumption of non-dissipative fireball 
dynamics. We find that F varies between about 200 and 
600 i^Vy^'^, while the three radii remain relatively stable 
with average values of rph = 2.5 ±0.8 x 10^^ C^^* Y^^^ cm, 
ro = 3.5 ± 1.7 X 10** y-'-^l'^ cm and = 1.3 ± 0.6 x 
10" ^-15/4 y-5/4 -yy^e especially note that the values of 
the photospheric and saturation radii are rather similar (see 
also Fig. |8}. Since the radiative efficiency of the thermal 
component scales as (rph/fg)"'^^'' (jMeszaros fc Re3l2000l ) 
this implies a high efficiency, as expected given that the 
photospheric component dominates the prompt emission. 

As discussed in section |4l the main uncertainty in these 
numbers is that we do not know how far the power law 
extends above the WAM energy band. In particular, if the 
power law extends to GeV energies the values quoted above 
for F and rph should be about two times higher, while the 
the values for ro and r^, will be smaller. In addition, all the 
above parameters have a weak dependence on the value of 
y, which is the ratio of the energy emitted in gamma rays 
to the total energy emitted by the burst. The value of Y can 
be measured through detailed aftergl ow observations an d is 
usually found to be less than 2-3 (e.g. lCenko et al1l2010l ). In 
the case of GRB 061007 it is hard to put a firm limit on Y as 
there is n o jet break seen in th e afterglow, but the preferred 
model of ISchadv erakl (|2007| ) gives a value of around 10. 
As in the case of a power law extending to high energies, 
this value would increase our estimates of for F and rph and 
decrease the values for ro and r^. 



5.3 The non-thermal component 

The power-law component in our model dominates at low 
and high energies and is consistent with having the same 
slope (s = —1.5) during the entire burst. Even though our 
energy range only extends up to 5 MeV it is interesting to 
compare this component with the high-energy power laws re- 
cently observed in some Fermi bursts. In particular, we note 
that a power law slope of — 1.5 has been observed in several 
cases l|Guiriec et al.ll2010a ') and that the isotropic energy re- 
lease of GRB 061007 is similar to the bright Fermi bursts in 
which this component has been observed. We also note that 
the power-law components that are required in addition to 
photospheric compo nents in BATSE G RBs have a preferred 
value close to -1.5 (iRvde et ahlbood ). 

There have been many suggestions for the origin of the 
extra component, including external sh ocks (|Ghisellini et al.l 
I2OIO!; iKu mar fc Barniol Duran' '2OIOI), hadronic processes 
(Asano et al. 2009; Razzaquc 20 1it), Compton u pscattering 
of a photospheric component (|Toma et al.ll2010l) as well as 
a com bination of different emission mechanisms l|Pe'er et al.l 
I2OIOI ). If the peak of the prompt emission is identified with 
a photospheric component it is also possible that the ad- 
ditional component is due to synchrotron emission arising 
in the region above the photosphere. The power-law slope 
of s = —1.5 supports this picture, as this is the slope ex- 
pected for fast-cooling synchrotron emission. However, in 
this case the power law break must be above the Suzaku en- 
ergy range, which implies an unreasonably large B-field (for 
a synchrotron peak at 1 GeV and F = 1000 the magnetic 
field is B ~ 10^'* G, where 7^ is the minimum Lorentz 
factor of the electrons). 
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Another option is that the power law is the inverse- 
Compton component of low-energy synchrotron emission. In 
this picture the synchrotron peak must be below the BAT 
energy range while the inverse Compton peak is above the 
WAM energy band. Such a large separation between the two 
peaks implies a very large efficiency for the electrons, which 
in turn rules out the standard internal-shock scenario. 

In order to evaluate models for the power-law compo- 
nent it is clearly crucial to determine their breaks, something 
which is now becoming possible with the Fermi satellite. So 
far, there is on e clear confirmation of a break in the case of 
GRB 090926A ijAckermann et al.ll201ll ). For GRB 061007 we 
expect a break of the power law in the LAT energy range 
(see section[4]), which should encourage the search for breaks 
in similarly bright bursts in the future. 
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